AEI-2009-111 
Imperial-TP-AR-2009-3 
UUITP-24/09 



From weak coupling to spinning strings 



Lisa FREYHULT a 0, Adam Rej 6 H and Stefan Zieme c H 

a Department of Physics and Astronomy, Uppsala University 
P.O. Box 803, S-75108, Uppsala, Sweden 

b Blackett Laboratory, Imperial College, London SW7 2AZ, U.K. 

Max-Planck-Institut fur Gravitationsphysik 
Albert- Einstein- Institut 
Am Muhlenberg 1, D- 14476 Potsdam, Germany 



Abstract 

We identify the gauge theory dual of a spinning string of minimal energy with spins S\ , S2 
on AdS§ and charge J on S 5 . For this purpose we focus on a certain set of local operators 
with two different types of covariant derivatives acting on complex scalar fields. We analyse 
the corresponding nested Bethe equations for the ground states in the limit of large spins. 
The auxiliary Bethe roots form certain string configurations in the complex plane, which 
enable us to derive integral equations for the leading and sub-leading contribution to the 
anomalous dimension. The results can be expressed through the observables of the sl(2) 
sub-sector, i.e. the cusp anomaly f(g) and the virtual scaling function Bi{g), rendering the 
strong-coupling analysis straightforward. Furthermore, we also study a particular sub-class 
of these operators specialising to a scaling limit with finite values of the second spin at weak 
and strong coupling. 
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1 Introduction and Summary 



Twist operators have so far played a major role in dynamical tests of the AdS/CFT corre- 
spondence in the planar limit. The main reason is their special scaling property at large values 
of the Lorentz spin. In particular, the anomalous dimension grows logarithmically with the 
spin, cf. pQ. This scaling behaviour is not a unique feature of the maximally supersymmet- 
ric Yang-Mills theory in four dimensions, but rather a special case of the so-called Sudakov 
scaling [2j common to many gauge theories. 

A typical representative of these operators in the 5 [(2) sector of N = 4 SYM is built from 
L complex scalar fields, Z, and M covariant light-cone derivatives, T>, acting on the scalar 
background fields within the trace 



The anomalous dimension of these operators occupy a band [3]. A distinguished sub-set 
among these is formed by the ground states, i.e. the lowest operators in the band, which 
enjoy several additional symmetry properties. For the two lowest possible values of the length, 
L = 2 and for the ground state of L = 3 , analytic expressions for the anomalous dimensions 
can be found at high orders in perturbation theory [U [5] . 

For twist-two operators these coincide with the maximal transcendental terms [6] of the 
known QCD results up to three-loop order, see [7] and references therein. At four-loop order 
the splitting functions of QCD are unknown, nevertheless the anomalous dimension of these 
operators may be determined (8J [9]. The result agrees with constraints from the BFKL 
equation [6]. Moreover, the result for M = 2 coincides with the explicit Feynman diagram 
computations of [10J. 

The ground states for L = 3 enjoy a similar solvability and the leading wrapping correction 
may be explicitly found [11J. In the special case of two excitations, the result is confirmed by 
the super-graph computation of [12J. 

For L > 3 it is unknown whether closed expressions for the anomalous dimensions can be 
found, even for the ground states. Nevertheless, the anomalous dimensions of the latter enjoy 
very interesting scaling properties in the limit M — > oo, 



The function f(g) is also referred to as the cusp anomalous dimension. It is conjectured to 
be independent of the length L and consequently not influenced by wrapping interactions. 
Thus, one can use the asymptotic Bethe equations to derive an integral equation [13], which 
allows to compute f(g) to arbitrary loop order. The scaling function resulting from the 
weak-coupling solution of this equation coincides up to four-loop order with the explicit 
perturbative computations of [H]. At strong coupling, the solution to this integral equation 
|15} IT6] leads to a remarkable agreement with the corresponding string theory results, see 
[TY[ [I8l [19]. Hence f(g) embodies the first known interpolating function of AdS/CFT. A 
special phenomenological interest in this object is due to its appearance in multi-loop gluon 
scattering amplitudes as well as in expectation values of certain Wilson lines. That is, the 
scaling function f(g) determines the leading 1/e 2 pole structure of the logarithm of gluon 
amplitudes [14J as well as the logarithmic growth of the anomalous dimension of light-like 
Wilson loops with a cusp [20] . as first noted in the strong coupling limit [21J. 

The virtual scaling function, B^g), appearing in (|1.2p explicitly depends on the twist L 
and it is less obvious that it remains unaffected by wrapping effects. The integral equation 



trV M Z L + ... . 




7sl(2) (L, M) = f(g)QogM + 7e - (L - 2) log 2) + B L (g) +0{^) . 



(1.2) 
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corresponding to Bi(g) has been derived in [22J, see also [23l [2¥J [25] . Interestingly, the solution 
to this equation may be related to the solution of the integral equation determining /(<?). 
This intertwines the strong coupling analysis of both functions and the methods developed 
for the cusp anomalous dimension may be directly applied also to the case of the virtual 
scaling function. The resulting strong-coupling expansion [22] is in perfect agreement with 
the string theory predictions at leading and next-to- leading order in A, see [26], suggesting 
that the wrapping interactions can be neglected also for the first finite-spin corrections. Also 
this quantity appears in the context of gluon amplitudes and Wilson loops. It enters the 
sub- leading 1/e poles as part of the collinear anomalous dimension, see |27) . 

In this paper we will go beyond the sl(2) sector. We introduce and investigate a gauge 
theory dual of a minimal energy spinning string configuration with two spins, Si and S2, 
in AdS§ and charge J in S 5 . The field content of these operators can be schematically 
represented by 

txV n+m i> m Z L . (1.3) 

The charges of the string are related to m and n through the identification 

Si = n + m— |, Si = m— |, J = L. (1.4) 

At weak coupling we extensively examine the limit m, n — > 00 with n/m = a fixed, in which 
S'l, S2 — > 00 with S1/S2 = 1 + a fixed. We start with the analysis of the corresponding nested 
one-loop Bethe equations. Surprisingly, the states with minimal length, L = 3, are again 
solvable and the respective Baxter functions may be found. While we extensively use this 
analytic solution for a numerical study of the behaviour of the one-loop Bethe roots in the 
large m, n limit, we defer its derivation to Appendix [A] We find that the auxiliary roots may 
be decoupled at the first two orders inm,n> 1, leaving a remainder in the main equation. 
Upon introducing the density of roots, this effective equation may be turned into a solvable 
integral equation. Subsequently, we use its leading one-loop solution as the starting point for 
the derivation of the all-loop integral equation for the density. The solution to this all-loop 
equation allows to determine the corresponding anomalous dimension. For the first two orders 
in m, n, we find 

7 (L, m, a) = I 7sl(2) (| L, m) + M l og (1(1 + a)(2 + a)) + o(^) • (1.5) 

This result is very surprising, since it suggests a deep relation between the one-spin and 
two-spin solutions at the first two leading orders. This unexpected link calls for further 
investigation on the string theory side. 

At leading order in spin the result ~ f/Q?) agrees with the energy scaling of spiky-strings 
|28[ [29] that consist of three arcs, each of them contributing a factor of h, see [30[ 131] . At 
finite order the results differ [261 E2] ■ 

Very recently it was shown by A. Tirziu and A. Tseytlin that in the case of a = the 
string theory calculation agrees with the formula (jl.5p at leading order in large A, see |33j. 
This gives first independent evidenc^] of absence of wrapping effects in our gauge theory 
calculation and justifies the use of the asymptotic Bethe equations. 

We furthermore examine the large n limit of fjl .3f) for arbitrary finite values of m > 2. 
The strong coupling solution cannot be expressed in terms of known observables. However, 
the same methodology as in [22] can be used, to solve the corresponding integral equation at 
strong coupling. At leading order in n this solution resembles the folded string configuration. 

1 It would be interesting to calculate the first quantum correction. 
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Figure 1: Dynkin diagrams of su(2,2|4) with different gradings [34] . 

2 Definitions and the string theory dual 

The highest-weight states corresponding to operators (II. 3p have the following Dynkin labels 
with respect to the upper Dynkin diagram in figure [H cf. [M] , 

(A , s\, s 2 , qi,p, q2, B, L) = (L + 2m + n — §, 2m + n — 1, n, 1, L — 3, 2, \,L) . (2.1) 

This suggests that the string theory dual of these operators is a spinning string with two 
spins S\,S% on AdS§ and charge L on the S 5 . Upon the usual 50(4) rotation, we can read 
of the spins directly from (|2.1|) 

S\ = m + n — i , S2 = m — \. (2.2) 

In view of the comparison with the string results, we will consider large values of m and n 
while fixing their ratio to ^ = a. For finite values of a both spins become large and 

| = l + «. (2.3) 

In particular, the case of a = corresponds to the symmetric case of the spinning string 
with equal spins. The well-studied case of the folded string with one large spin, on the other 
hand, corresponds to a = 00. In this publication, we will only study the gauge theory states 
with minimal anomalous dimension, which correspond to minimal energy states on the string 
theory side. 

In contradistinction to the sl(2) sub-sector, the Bethe equations corresponding to the op- 
erators (|1.3|) are nested. The number of nesting levels depends on the choice of representation. 
The minimal number of levels is equal to three. However, for the sake of convenience, we 
will mostly use a non-minimal representation corresponding to the lower Dynkin diagram in 
figure [TJ 

A subset of (|1.3|) with L = 3 has already been studied in the literature, cf. |35j. It was 
found experimentally that the corresponding one-loop anomalous dimension is given by the 
following closed formula 

7 n , m = 4ffi(s_i) + 4fr 1 (m + a)+4 J ffi(f + 5)_4fT 1 (-i) . (2.4) 

Here, Hi(N) denotes iV-th harmonic number. We will prove this formula in Appendix 1X1 
by constructing an explicit one-loop solution for L = 3. This is a counterpart of the L = 2 
solution found in the s[(2) sector [36]. The corresponding Baxter functions (see (jA.28|) . ()A.29|) 
and (|A.34j) ) become quite complicated and are not hypergeometric orthogonal polynomials 
anymore. Nevertheless, the explicit one-loop solution enormously facilitates numerical studies 
of the Bethe solutions at large values of m and n. 
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3 Large spins solution at one-loop order 



In this section we will analyse the operators (jl.3p in the limit m, n — > oo and n/m = a = fixed 
at leading order in perturbation theory. As in the case of the s[(2) operators [36J, the leading 
solution does not depend on the length L. The minimal set of equations at one- loop consists 
of two coupled nesting levels 

(. i \ L n+2m . m . % 

u A,k + 2 j _ TT u 4,,k ~ Ujj ~ I -j-r U^ k - U 5 J + g 

m . n+2m . j 

l = TT U5 ' k ~ U5J ~ % TT u ^ k ~ UA ' j 5 (32) 

From the numerical studies of the analytic solution for L = 3 presented in Appendix |A] we 
infer that the roots 115 form strings 

U5 — ¥ ±i/2, ±3«'/2, ± . . . m even, 

U5 — >• 0, dbi, ±2i, ± . . . m odd. (3-3) 

However, it is incorrect to assume, even at the leading order, that all u§ will exhibit this 
behaviour. We thus introduce an effective cut-off c(a) such that 

U5 — >• ±i/2, ±3i/2, . . . , ±i(m/c(a) — 1/2) m even, 
u 5 -»0,±*,±2i,...,±i(m/c(a)-l/2) m odd. (3.4) 

In the following we will explicitly fix c(a) in the limits m,n — >■ 00. We have also checked 
numerically that the remaining roots, i.e. not belonging to the effective strings (|3.4p . scale 
as ~ m? at large values of m and thus are not relevant for the leading and the sub-leading 
order. In the large n limit and for finite values of m the strings (|3.3p . as we will show in what 
follows, become exact and c(oo) = 2 as expected. In either case, one can completely decouple 
the W5 roots. The effective equation for the middle-node roots thus takes the following form 

(j \ L / 1 im \ m(a+2) 

«4,k - 1 / \U 4 ,k 7 = i U 4,fc - U4J + » ' 

In the large m and/or n limit this equation may be turned into an integral equation along 
the lines presented in |36j . Explicitly, one obtains 

arctan(2«fc) = — r^( u ) — 2 f du'p(u') arctan(u — v!) 

m[a + 2) Ju a) 



m(a + 2) 



Here, n(u) is the function for the mode numbers. It is related to the density po(u) through 
n(u) = -— sgn(u) - i + / du'p(u') . (3.7) 

Z 1 J-b(a) 
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3.1 The a = solution 

Upon rescaling the momentum-carrying roots = 2m and the corresponding density 
Po(u) = 2mpo(u), the leading part of the equation (|3,6p is given by 

/• 5(0) dn(u') 

37rsgn(u)-2 arctan(2cM) + 2-/ du' _ K _,. = 0. (3.8) 

J-b(o) ( u ~ u ) 

The solution to this equation is given by 

/ /i " 2 i V 45 (o) 2c2 -^ 

6(0) 2 2fe(0)c 



^ ) = ^" K (2^))-27 l0g 



. V45(0)2 C 2 + 1 _ / _ _^ 

\ 2b(0)c V 6(0) 



(3.9) 



Here, with p K ( u ) we have denoted the one- loop density of the Bethe roots corresponding to 
the ground state of twist-L operators in sl(2). The normalisation condition 

r5(0) 

1 = / dup (u) , (3.10) 

J-b(0) 

relates the boundary of the root distribution to the dimension of the effective strings c. 
Explicitly, one finds 

6(0) = — (\/9-4c + 4c 2 + 6c-3). (3.11) 
16 c 

To determine the constant c we proceed as follows. The density po(u) has the following large 
m expansion 



P ^ U) = ^{n l ° gm + C ~^ l0g(n V + ' ' ' ' (3-12) 
with the constant term C given by 

C = ^ log (4 6(0)) - ^ log ^2 6(0) c + ^4 6(0) 2 c 2 + l^) . (3.13) 



One may thus split the density in fj3.6j) as follows 

p(u) =p («) + r(n), (3.14) 

and use the expansion (|3.12p to obtain a leading integral equation for r(u). Upon Fourier 
transformation one finds 

f(t) = — | 3 € m M - ^ ^— | . (3.15) 

w 2m 1 1 — e— 1*1 \t\ l + e-l*l/ 2 J V 1 

One can now use (j3.14p to calculate the one-loop anomalous dimension. After straightforward 
integration one finds 

70 = 121ogm + 127 E - 8(L - 3) log 2 + 6 vr C. (3.16) 

This should be compared with the large m expansion of the L = 3 analytic result (12.41) 

7o i=3 = 121ogm + 12 7E . (3.17) 

Thus, the constant c is determined by the condition C = in conjunction with (|3.1ip . 
Numerically one may determine c = 2.83181(...), which means that approximately jq of the 
U5 roots form effective strings. We have checked numerically that 6(0) obtained by inserting 
this value of c is consistent with the scaling properties of the largest 114 roots. 
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3.2 The solution for general a > 

The limit m — > oo, n — > oo with ^ = a = fixed is a simple generalisation of the a = case 
discussed above. The roots U5 form again strings, and one expects the cut-off parameter to 
depend on a, i.e. c = c(a). Upon rescaling the roots by n + 2m = m(a + 2), one derives the 
following integral equation 

/•&(«) p n (u' a ) 

3vrsgn(u) -2 arctan((2 + a) c(a)u) + 2 4 du' ,_ K ' ; = 0. (3.18) 

y-S(a) («-0 

The solution to this equation is given by 

/ R ^ 1 Vfc(») 2 («+2) 2 c(a) 2 +l ' 
V 1 fe(a) 2 b(a)(a+2)c(a) 



yjb{a) 2 (a+2) 2 c(a) 2 +l 



P " 0( "' a) = 2^ K (2^))-27 log j , _ 

\ b(aj(a+2)c(a) \J 1 b(a 

The normalisation condition yields the relation 
1 



(3.19) 



= ^ \ — 7 — \ ( -3 + 3 (2 + a) c(a) + ^9 - 2 (2 + a) c(a) + (2 + a) 2 c(a) 2 ) . (3.20) 

o(2 + a) c(a) \ J 

The constant c(a) may be determined by the same procedure as in the a = case. The 
denisty po(u) has the following large m expansion 

3/2 2 

Po(u,a) = — r -logm + (7(a) , log(n 2 ) 

2m(2 + a) \tt 7r z 

with 

C(a) = - log (2 5(a)(a + 2)) - A log f 6(a)(a + 2)c(a) + \/b(a) 2 {a + 2) 2 c(a) 2 + . 

7T 07T V / 

(3.21) 

Splitting again the leading density as 

p(u, a) = p (u, a) + r(u, a) + . . . , (3.22) 

one determines f(t, a) to be 

r(t, a) = — — ^ 3— nr - A • (3-23) 

v ' ; m(2 + a) ^ l-eH*l |t| 1 + e -l*l/ 2 y V ; 

This immediately leads to 

70 = 12 log m + 12 7e - 8(L - 3) log 2 + 6 vr (7(a). (3.24) 
Expanding (12.4j) for general a, one finds 

70 = 121ogm + 12 7E +41og(i(l + a)(2 + a)) . (3.25) 
One thus concludes that 

3 vr C(a) = 2 log (|(1 + a)(2 + a)) . (3.26) 

The above formula in conjunction with (I3.2ip and (I3.20p determines c(a) uniquely. The 
behaviour of c(q) as a function of a is shown in figure [2j Clearly, the effective strings become 
exact in the large a limit 

lim c(a) -> 2 . (3.27) 
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2.6 

c(a) 24 



a 

Figure 2: The plot of c(a) as function of a. 



3.3 The a = oo solution 

In view of (|3.27|) , one concludes that for a — > oo all roots u§ form strings as in (|3.3p . Rescaling 
the main roots as U4 = n U4 and taking the large n limit we find 

■b/n ( u ~ u ) 

The above equation, together with the momentum constraint which fixes b to n/2, is solved 
by the known sl(2) density, p K {u), found in [37]. This formally corresponds to taking a — > 00, 
although the latter is not unique due to the fact that any dependence on finite values of the 
spin 52 is lost. 



4 The all-loop equations 

Beyond the one-loop order the spectral equations for the operators f) 1 .31) can be derived from 
the full system of the asymptotic Bethe equations conjectured in [34] . In what follows, we 
will work with a set of four coupled Bethe equations 

/ + \ L n+2m - _ + i_„2/ + - 

-=H =1 —L ~ 9 ' + (T«4,*»H: 

I ry I -A. -M- ry 1 ry I ft £ I ry ry ' 

\ X 4.fc / s-/.u x 4.fe x 4,? 1 iJ I^ai-^aa 



x TT TT 

m— 1 . /0 n+2m „— 

1 771 / TTL — 2 / 

_ TT U<5,k ~ U(5,j + j TT ^6,fc ~ ^5,j ~ V2 yr Ugj ~ UjJ ~ E/2 

M w 6)fc - u 6 j -ill n 6jfc - it 5 j + i/2 11 n 6)fc - it 7 j + i/2 

m— 1 ■ / n n+2m i 9 /„ _— 

= TT ^7,fc ~ «6j - »/2 TT 1 ~ g /X 7 ,fc^ 4 ,j f , 

1 = i 1 «7,fc " U 6>j + i/2 11 1 - 5 2/x 7 , fc 4, ' 



7 



Our convention for the coupling constant is g 2 = g 2 M N c /(16ir 2 ). The deformation of the 
spectral parameter reads x(u) = \{u + y u 2 — 4<? 2 ), with the conventional notation x ± = 
x(u± |). The form of the dressing factor a 2 is given in [13] . Note that at one-loop order, the 
generating functions of the uq and uj roots can be immediately obtained from the U5 roots. 
Consequently, these roots inherit the behaviour of the M5 roots and also form effective strings. 
We will denote the effective cut-off parameter for the uj roots by d(a). 

In the large m, n limit, the leading positions of the inner 115 roots are again given by ()3.4p . 
The quantum corrections to these roots vanish as m,n —¥ 00. The outer roots, on the other 
hand, grow very rapidly. The same is true for the auxiliary roots u-j. Thus, effectively, the 
system (|4.ip reduces to 



x k 



n+2m 



9 2 lx\x- 



II Uk _ Uj + iy 1 _ g 2/ Xk x + j 

From this set of effective Bethe equations we will derive an integral for the fluctuation density, 
a(g,u,a), 

p(g,u,a) = p(u,a) — —a(g,u,a). (4.3) 

m[a + 2) 

In what follows, we will treat the cases a = 0, a > and a = 00 separately, even though the 
first two directly interpolate between each other. The case of a = 00 requires a more careful 
analysis. 



4.1 The equal spin case a = 

Knowing ( I 3.12p and (|3.15p , one can ea 
transforms of the fluctuation density a(g, t, 0) 



Knowing (|3.12p and (|3.15p , one can easily derive the integral equation for the Fourier-Laplace 



t 

e~- 1 



L 



&(g, t, 0) = -j—. r (§ log m + §7 E - - 3) log 2) if (2^, 0) - ^ (J (2gt) - 1) 



1 r 00 / 3 r — 3 \ 

+ 2 / " ^Tl) (K(2 9 t,2gt')-K(2gt,0)) 

-Ag 2 / dt'K(2gt,2gt')a(g,t',0) . (4.4) 

It is written in terms of the usual integral kernel K(t,t') = Ko(t,t') + Kx(t,t') + Kd(t,t'), 
with the parity even and odd components given respectively by |13j 



00 



n=l 

rr (+ +'\ t'Ji(t)J (t') -tJ (t)Ji(f) 2 ^ 

n=l 



2 We define the Fourier-Laplace transform of the fluctuation density by a(t) = e 2 due ltu a(u). 
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The dressing kernel Kd(t,t') is a convolution of the even and the odd part 



f°° t 
K d (t,t) = 8g 2 dt?K 1 (t,2glf , )- p} -—Ko(2ge , ,l/). (4.6) 
Jo e* - 1 

The anomalous dimension corresponds to the value of a(g, t, 0) at the origin 

7(L,m) = 16 5 2 <r( 5 ,0,0). (4.7) 

By comparing (14. 4p with the corresponding equation for the sl(2) sector |22j . one easily infers 
that terms proportional to K(2gt, 0) in the first line of (|4.4|) give rise to the cusp anomalous 
dimension f(g), while the remaining terms yield the virtual scaling function -£>2(fi0 obtained 
in [22], and the first generalised scaling function e\(g) of [381 US]- Hence, the anomalous 
dimension is given by 

7 (L, m) = \f{g) (logm + 7e ) + (£ - 3)ei( 5 ) + ^ 2 ( 5 ) + o( ' 



2 J wy v ° ' ty v ' w/ 2 w Vlogm 

= | WSA ™)+o(i^). M) 

Remarkably, the large m anomalous dimension is up to the order 0( lo * ) proportional to 
the anomalous dimension of the twist operators. This can be traced back to the "screening 
properties" of the us roots, cf. (13. 4p . It should be noted that although we have assumed 
a = 0, one does not need to put n = 0. In contrary, (|4.8p is valid also in the latter case 
and the first finite n corrections are sub-leading. We expect, in similarity to the case of sl(2) 
operators, that the wrapping corrections will not affect the first two orders in the large spin, 
S = 2S±, expansion. The findings of [33J confirm this hypothesis. 

4.2 The case of a > 

It is straightforward to repeat the computation of the preceding paragraph for a > 0, the only 
difference being an additional one- loop term in the energy, cf. (|3.25|) . The resulting equation 
for a(g, t, a) takes the following form 

a(g, t, a) = [ (| logm + | 7e - (L - 3) log2 + \ log (|(1 + a) (2 + a))) K(2gt, 0) 

T 1 / 3 T — 3 \ 

-_(j o (20t)-l) + - ^ ^(_-^_)(jr(2 l7 t 1 2^)-^(2^0)) 



oo 







4 5 2 / dt'K(2gt,2gt')a(g,t',a) . (4.9) 



The anomalous dimension can be easily found 



1l(9, m) = l %m {lL, m) + ^ log (1(1 + a)(2 + a)) + O (^ ) . (4.10) 



It is noteworthy that the dependence on a is logarithmic and that the corresponding prefactor 
is again proportional to f(g). It would be interesting to understand from a string theory 
perspective why the energy of the general two-spin solution so closely resembles the one-spin 
solution. 
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4.3 One large spin limit 

In this section we will discuss the case of n — > oo and finite values of m. It turns out that, 
contrary to the a = case, the sub-leading correction exhibits an interesting dependence on 
m. 

The one- loop leading solution has been discussed in section 13.31 and is fully equivalent to 
the one-loop problem for the ground states in the sl(2) sub-sector, cf. [HB]. Proceeding in the 
same spirit as in the previous sections, we split off the leading density 

p{u) = p (u) +r(u) (4.11) 

and derive a leading equation for r(u), which subsequently may be solved by Fourier trans- 
formation 

Using this expression, it is straightforward to calculate the one-loop anomalous dimension for 
the first two orders in n 

E = 41ogn + 6 7e + 2M m j 1 )-<L- 3) log 2 . (4.13) 

The derivation of higher-loop corrections goes along similar lines as in the preceding sections. 
Upon defining the fluctuation density by 

p(g,u,m) = p(u,m) a(g,u,m) , (4.14) 

n 

one derives the following closed integral equation for m>2 

*(g, t, m) = -J-^ [ (logn + | 7e + ~(L-3) log 2) K(2gt, 0) - ^ (J (2gt) - 1) 

i poo . o -(m+l)t'/2 T _ o . 

+ 2 J, dt '{^T - i _ e -t> ~ TT^Tl) (K(2gt,2gt')-K(2gt,0)) 

of 00 1 f°° ,, m-l 

-Ag 2 dt'K(2gt,2gt')a(g,t')-- R^gt^gi'^ 2 . (4.15) 
Jo 1 Jo 

The resulting anomalous dimension scales logarithmically with n 

7(L, n, m) = I6g 2 a(g, 0, m) = f(g) log n + . . . , (4.16) 

while the finite-spin corrections depend explicitly on m, as can be directly inferred from (14.151) . 
At weak-coupling, one can easily determine the perturbative expansion at the first few orders 



7 (L, n, m) = f(g)(logn + § 7e + ^o(^) - (L - 3) log 2) - 2g 4 (M m ^ 1 ) + 2(21 - 4L)C(3)) 
+ 1/ ^ 4 (-±i) + 2^(2^1) + 24Vi(=^)C(3) + (^tfC(3) 
+8(6 - L)vr 2 C(3) + 72(31 - 7L)((5)) + ... , (4.17) 



For the choice of parameters m = 2 and L = 3 we find perfect agreement with the large 
n expansion of the anomalous dimension up to four-loop order which is available in [10] . 
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However, since the g 8 contribution to ()4.17j) is quite lengthy, we merely give the first three 
loop orders. 

Interestingly, equation (I4.15D can also be solved at large values of the coupling by making 
use of the strong coupling expansion for twist operators in the s[(2) sector [15]. We defer the 
strong-coupling analysis to Appendix [B] and only present the final result 

7 {L, n,m) = Ug ^_ i og _ + 6 </ (log2 — 1) + (1 — L) + + - 

V 7r / g m — 12 



5 Outlook 

The appearance of the cusp anomalous dimension f{g) and the virtual scaling function B^{g) 
beyond the st(2) sector in (jl.5p is quite remarkable and certainly promotes their universality. 
For further sub-leading corrections in spin, the endpoints of the effective condensate may not 
be enough to completely determine the corresponding contribution. In this case the remaining 
roots should also be taken into account and it is questionable if the scaling still resembles the 
behaviour of twist operators. 

It will be quite interesting to investigate if the generalised scaling function f(g,j) of the 
sl(2) sector |23j also appears in the refined limit S — > oo, L — > oo with j = j^|hj fixed. 

In similarity to the known solvable cases of twist-two and three operators, it would be 
very interesting to see whether it is possible to construct higher loop contributions to the 
anomalous dimension (|2.4p for general values of m and n. A first step in this direction has 
been made in [3U] , where the case of m = 2 was analysed. 

Furthermore, the decoupling procedure described in detail in Appendix [A] is based on 
iteratively splitting one bosonic node of the corresponding Dynkin diagram into two fermionic 
ones. Although this is straightforward at the level of the equations, it remains obscure to us 
what the corresponding algebraic interpretation might be. 
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A The analytic one-loop solution 



In this section we will study the following class of length-three operators 

t?V n+m V m Z 3 . (A.l) 

For m = this set reduces to twist-three operators of the sl(2) sub-sector. Also the case 
of m = 1 is redundant as (lA.ip then corresponds to descendents of twist-two operators. 
Therefore in what follows we will assume m > 1. A subgroup of these operators for m = 2 
has been studied in [ID] up to four-loop order. At one-loop order, on the other hand, a closed 
expression for the anomalous dimension of the ground states has been conjectured in |35j for 
any value of m and n, see (|2.4p . It is rather straightforward to prove this formulsH using the 
analytic solution provided below. 

The excitation pattern for the higher-loop Dynkin diagram in the upper part of figure Q] 
reads 

(K 1 ,K 2 ,K 3 ,K 4 ,K 5 ,K 6 ,K 7 ) = (0, 0, n + 2m - 1, n + 2m, n + 2{m - 1), m - 1, 0) , (A. 2) 

where K u denotes the excitation number of the v-th node of the Dynkin diagram. After a 
dualisation of the 113 roots the corresponding one-loop system of equations describing this 
class of operators is given by 

(, i \ 3 n+2m-2 j 

m—1 j n+2m j 

1 _ TT ^5,fc ~ M6,j + 2 TT n 5,fc ~ ^4,j ~ g ^ 4 ) 

J= l U5,k ~ U 6 j - | U 5ik - U 4 ,j + § 

m—1 . n+2m— 2 j 

This set of Bethe equations is valid for m > 1. In the following section we will solve it exactly 
thanks to a hidden recurrence relation between roots for different value of m. 

A.l The one-loop recurrence 

We start by dualizing equation (jA.3p . The system of equations ()A.3[) - ()A.5[) reduces to 

1 -\ 3 n+2m— 2 . m—1 j 

m—1 . ra+2m— 2 . j 

3 M. Beccaria, private communication. 
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The derivation of the hidden recurrence is based on the observation that the last equation is 
equivalent to a system of two coupled equations 

m— 2 i n+2m— 2 . j 

1= JJ U 6,fc ~ U7,j ~ 2 j-J u 6,fc ~ U 5J + | ^ ^ 



m—1 



2 



J=l n 7,fc - U6,j ~ 2 

where we have introduced a new set of auxiliary roots uj. We now introduce the Baxter 
function for the uq roots and their dual counterpart uq 

m— 2 n+2m— 2 m— 2 ra+2m— 2 

R( u )=Yl(u-u 7d -^) jj ( u - u 5,j + ^) - Yl(u-u 7J + ^) jj («-«5j-5) 
j=i j=i j=i j=i 

m— 1 n+2m- 4 

= c 6 (u - u 6 j) JJ ( u ~ S 6,i) • (A. 10) 

It is straightforward to derive the following two relations 

r>t i i\ n+2m— 2 . m—1 . j n+2m— 4 ~ . j 

JH^5,fc + gj _ TT ^5,A: ~ ^5J + 2 _ TT ^5,fc ~ ^6,j + g TT M 5,fc ~ Uqj + g 

#(^5,fc - §) J= i u 5,fc - U5,j - i jJi U^ k - Uqj - | j-Ji U^ k - Uej - I ' 

t-,1 i i\ m—2 . m—1 i n+2m— 4 ~ . j 

Jt(,^7,fc + gj _ TT u 7,k ~ u 7,j + 1 _ TT n 7,fc ~ ^6,j + g TT u 7,fc ~ ^6,j + g 12 ) 

With the help of (|A. 1 1[) we now rewrite (|A.6|) as 

. • \ 3 n+2m-4 ~ j 

^5,fc + M = TT u 5,k ~ U 6J - | 

U5,k-i) fJi u 5ifc -n 6ii + ^' 
and once again dualize this set of equations by defining the polynomial Qs(u) as 

n+2m-4 n+2m-4 n+2m— 2 

Q 5 (u) = (u + i) 3 Yl ( n ~ + §) - (« - i) 3 JJ - «6,j - I) = c 5 JJ (u-u 5 j). 

j=l j=l j=l 

(A.14) 

The function Q§{u) obeys the relation 

t ~ j \ / ~ , a . \ 3 n+2m— 4 ~ ~ , . n+2m-2 ~ , 

Q5{Ue,k + |) _ / ~»6,fc + 2 Z I TT U G,k - U6,j + \ _ TT ^6,fc ~ U 5 J + | ^ 
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Since the uq roots also solve (!A.8j) . we can decouple the M5 roots. Likewise, we use (1A.12|) to 
rewrite ()A.9j) in terms of 5g . The resulting set of equations thus reads 

(~ 1 3 A 3 n+2m— 4 ~ ~ . m-2 ~ , i 

m— 2 . n+2m-4 — j 

fj[ u 7,k ~ u 7>j + i JL1 U7 fe _ ^ _ i • 

Comparing the two systems of equations, namely (|A.6|) - (|A.7|) with (|A.16P - (|A.17p . we note 
that the value of to has been lowered by one, while the spin of the representation increased 
by \. Clearly, this procedure can be applied recursively until all second-level roots vanish 
and get absorbed into the spin representation of the first-level roots. Thus, after (to — 1) steps 
one can decouple ([A.7|) from (|A,6j) and one is left with a single system of equation^. Before 
using this recursion to solve the system (|A.6|) - (|A.7|) . we will investigate the Baxter functions 
appearing in the intermediate steps. 

Suppose that one has repeated the aforementioned procedure £ times. The intermediate 
equations then read 

\U k {l-\- 2 )l/ j=1 U k Uj +1 j=1 U k V- 2 

1= TT v Ji V J ! TT A u i + 2 (A19) 



j=i u fc ~ v j + l j=i v k ~ u j 

j+k 



with the initial values 



u[ = U5,k , = u 6jk , (A.20) 

(°) — C 1 ) — C A 01 A 

«fc = «6,fe , v k = U T,k • ( A - 21 ) 

A single step of the iteration relates the polynomials 

n+2{m-l-t) 

P £ (u)^ [] (u-uf), (A.22) 
3=1 

with consecutive values of I through 

P e (u) = ( u + (l + l)i) 3 P e+1 ( u +i)-( u -( 1 + l)i) 3 Pe+1 ( u - i) . (A.23) 

The initial polynomials (jA.lOj) and (|A.14p are respectively the second and the third member 
of this family, i.e. Qe(u) = Pi(u) and Qs(u) = Po(u). It should be clear that the first element 
is given by 

n+2ra 

P-i(u) = Qa(u) = c 4 J] (u - U4j) . (A.24) 



4 Please note, however, that in the penultimate step one should not perform the splitting, but rather dualize 
the nested set of equations directly. 
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A general solution to the recurrence relation (1A.23|1 is given by 

w -D-<)n(--^)'n(- + ^)V(. + (=-»). 



k=0 x y 3=1 v ' j=l 



(A.25) 



with n being an arbitrary positive integer. 

As already mentioned before, for I = m — 1 equations (|A.18[) - ()A.19|) decouple and one is 
left with a single equation 

/ (m-l) + n^A ^ n M (m-1) _ u (m-l) _ ■ 

1 (m-l) m +l • J ~ 11 (m-l) (m-l) . - ' (A.26) 
\ u k ~^t % ) j=l u k ~ u 6 + l 

It is noteworthy that these are the Bethe equations of a non-compact sl(2) magnet in the 
spin-(— ^j^) representation. Thus, we have proven an equivalence between both systems 
noticed in [35]. The corresponding Baxter equation can be solved exactly for the ground 
state. Please refer to App endix I A . 2 1 for further details. The solution is the Wilson polynomial 



F n ,m(v) =P m _i(w) = c m _i Y\(u ~ u 



(m— 1) \ 
'3 

3=1 



: 4-^3 



' 2 ' 2 ^ 2 ' 2 2 * U 

J. -|- 2 , i -|- 2 , J- i 2 



(A.27) 



Plugging this into ()A.25j) . one finds an explicit solution for the primary roots 

m , \ k n m—k „ 

Q 4 (n) = ^(-l) fc r)n(--^) n(- +M ^) ^»,m («+(?-*)<), (A-28) 

fc=0 ^ ' j=l 3=1 

while the U5 roots are generated by 

m— 1 / \ fc m— 1— ft 

Q5(«)=^(-l) fe ( m r" )n^-^) 3 II (^+Ji) 3 i ? n,m(n+(^ 1 -A ; )i) . (A.29) 
fc=0 ^ ' 3=1 3=1 

The Baxter function for the «6 roots is a part of a different recursive scheme. Defining 

m-l-£ 

K t {u)= [J (u-vf), (A.30) 
j'=i 

we find that Kg obeys the following functional relation 

K t (u) = f e (u) K e+1 (u + |) - /,(«) tf m , (A.31) 

with 

AM = -^ff and £(«) (A.32) 

"<+iM %i(«) 
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p = -l 



«4 



"5 



U 6 

■o 



U5 



>=°;o— o 
»= +i o— o: 



uCp) ' v (p) 

o— o; 



\ u (pi i; \ t .(pi i; 

» +i O— O 



w 5 

D- 



"6 



«7 



;o— 



Figure 3: A schematic representation of the decoupling procedure. Dotted lines indicate 
the bosonic splitting and the dualization of the interjacent roots. The curved lines denote 
two kinds of recurrence relations P and K, respectively. The exact one-loop solutions Q4 
and Q5 belong to P iteration, while Qq is part of the K recurrence. 



A general solution to equation (IA,31j) is presented in Appendix IA.31 We fix the boundary 
conditions by specifying 

m— 1 

K m _i = 1 and K (u) = Qe(u) = JJ (u - u 6 j) . (A.33) 

i=l 

Thus, we find that Qq is given by the following expression 

m-2 



'try / k \ 
Qe{u)= [[ fk ( u+ 2 

1 — n V / 



j s - 2(r - s) 



jr-l 



k=0 



k=0 

m—1 m—2 jr-1—1 r , 

+E(- 1 )"EE- £ n&(-+ 

r=l ii=0 j2=0 >=0 s=l v 

>nn'4" ^- 2| 7 8+1) -) n/.(."i 

s=2 fc=j 3 +l ^ ' fc=ii+l ^ 

with f p (u) and / p (it) given in ()A.32p and 

/ m — 1 



n fk « + 



— ?. 
2 



(A.34) 



m—1- 



p e ( U )= (-1) 



fc=0 



A- 



X-F n ,m ( u + ( r "^— 1 — - - k I i I , 



(A.35) 
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as follows from (jA.25j) with £ + n = m — 1 . 

The complete solution of the one- loop problem is hence given by (|A,28|) , (|A.29P and (|A,34p . 

A. 2 Solution to the higher spin magnet 

The Bethe equations in (IA.26h correspond to an integrable non-compact XXX m+i magnet. 

2 

The Baxter equation associated with XXX^ S magnets takes the following generic form 

(u + is) L Q(u + i) + (u- is) L Q(u - i) = t L (u)Q(u) . (A.36) 

For the special case of L = 3, the transfer matrix of the ground states can be determined 
explicitly 

t 3 (u) = 2u 3 +q 2 u + q 3 , q 2 = -(n 2 - n + 6 n s + 6 s 2 ) , q 3 = . (A.37) 
The solution to (lA.36j) with (lA.37j) can be found by noting that the Wilson polynomials Wk , 
Wk(u 2 ; a, b, c, d) _ ( —k, k + a + b + c + d— 1, a + iu, a — iu 



1 , (A.38) 



(a + b)k (a + c)k (a + d) k — \ a + b, a + c, a + d 

satisfy the following difference equation (see e.g. |42j) 

k(k + a + b + c + d- l)y{u) = B(u)y(u + i) - [B{u) + D(u)]y(u) + D(u)y(u - i) . (A.39) 
Here, y(u) = Wk(u 2 ; a, b, c, d) and 

, (a — iu)(b — iu)(c — iu)(d — iu) . (a + iu)(b + iu)(c + iu)(d + iu) 

n(U) = ; : , V(U) = ; r . A.4U 

v ; 2iu(2iu - 1) ' v ' 2iu{2iu + 1) v 1 

It is easy to check that (lA~36l) is embedded in (lA~39l) upon the identification 

a = — , b = c = d = s and k = — . (A. 41) 

2 ' 2 y J 



Hence, the solution is given by 



(A.42) 



For the special value of s = one finds (|A.27h . 

A. 3 Solution of the recurrence 

The class of functional equations 



M u ) = f P ( u ) A v+i ( u + - f P (v) a p+i - 1 J ( A - 43 ) 
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is solved by 

n-l 



A p (u) = Y[ fp+k (u + - ij A p+n (u+ 7 ^ 

k=0 ^ ' 

+E(-D r EE-E nw ("+ ■" 2 ' ') n /**(«+§') 

r=l j 1= j 2 =0 >=0 s=l ^ ' fc=0 ^ ' 

n'n" 1 w(~+ ' i " 2(r 2 " s+1) i) ff 

o— 9 i A J_1 \ / I i_ J_1 \ / 



r=l j 1= j 2 =0 >=0 s=l v ' fc=0 

r 1 / , „/ , -i n n n—1 

X 

s=2 k=j s +l ' " ' fc=ji+l 

/ n — 2 r \ 

xA p+ „ u+ — - — »J . (A.44) 
The proof is by induction. 

B The strong coupling limit 

In order to analyse the sub- leading contribution to (|4,16|) we decompose the density a(t) into 
its parity even and odd parts 

Since the kernels Kq and K\ in fj4.5f) are given by the sum over Bessel functions, the functions 
7± take the form of a Neumann series 

oo oo 

7+ (2 5 t) = 2 ^ 2nJ 2n {2gt) l2 n , 7-(2<7*) = 2 ^(2n - 1) J 2n _i(2 5 t) 72n _ 1 . (B.2) 

n=l n=l 

Using this decomposition of the fluctuation density it is possible to rewrite the integral equa- 
tion (j4.15|) as an infinite system of equations with n > 1 

dt f 7+(t) 7-(*) > \ , m £ , , 1 r, j2n(2gt) _ t m=i 

J2n\t) = o 1 «2n — 7: / «t 77— 6 2 , 



o t \i- e -t/2g e t/2g-iJ znx ' 8ng in 2 J 2gt 



oo 



/o t \ 1 - e-*/ 2 f e*/23 _ i 
where the term /i n = h n (g) is given by the expression 

Since the left hand side of the equations (1B.3|) is the same as in the case of the BES equation, 
we expect to be able to express the solution in terms of the solution to the BES equation. 
For this purpose we introduce a new parameter j, which interpolates between the system 
corresponding to the BES equation and (|B.3[) 

dt ( 7+ (t,j) 7 _(t,j) ^ jL j f°° ^ J 2n (2gt) _ t ^i 

J2n[t) = \-Jtl2n-- / dt — e 2 , 







'0 



t \l- e -t/2g e t/2g_ 1 J ^> 8ng <•> 2 J 2gt 

+ ^7^ rH2™-iW = ^2n-i + -(l-j)^,i. (B.5) 



i \ l_ e -t/2 fl e t/2g_i ^-x\ , j zn-x ^ 
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Setting j = gives back the BES equation while j = 1 corresponds to (|B.3j) . Multiplying 
the first equation by (2n)7 2n (i, j') and the second by (2n — l)7 2n -i (t,j'), summing over all 
n and finally subtracting the two equations leads to a left hand side that is symmetric under 
exchange of j and j', see [38] for details. Using this fact and setting j = and j' = 1, we find 

1l ( Sl .i4r*(=--' i -"" , )(7#^+ t+m) 7i(s - o) 



4 Jo V ; V(e' /2;? - (e-*/ 2 a - l)gt {eWs -l)g 

L-3 [°° , ( 7-M) , 7+M) 7iM) \ mfil 

4 7o V(e*/ 49 + l)^ (e-'^ + l)^ (e*/4fl + ■ 

The finite order correction is then given by 16 g 2 71 (g, 1). A change of variables as in |15| . 

27±(*,0) = (l-sech(^))r ± (t,0)±tanh(^)r T (t,0), (B.7) 



leads to 

1 , , :; 

16g 



7i(5, 1) = 7^2 ( L " 3 ) £ i(3) + 7i(5, 0)(L - 3) log 2 + -5 2 ( 5 ) 



4jr- dte ^^(^cr + (*,o) + r^o) ) + ^^) . (b.8) 

At this stage, we make use of the solution of the BES equation obtained in [TS] , 

00 00 

r+(t,o) = ^(-i) fc+1 J 2fc (t)r 2fc , r_(t,o) = ^(-i) fc+lj 2 fe -i(*)r 2fc _i , 

fc=0 A;=0 

where the coefficients Tk are given by 

1 00 1 

^-4 0) + E^ferr+«), (B.9) 



7-' : 



1? 



{p) _ T(m + p-±) w _ (-i)Pr(m-i) 

2m r(m + i)r(i)' ^ r(m + i- P )r(i)" 1 j 

The prefactors Cp explicitly depend on g and can be determined from the so-called all-loop 
quantization condition of |15j . Bearing that in mind, we find for the integral in (IB.8j) 




CB.11) 
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Subsequently, using that Tq = Ag^i(g,0) and T_i = 1, we can recast 1(g) as 



oo / ■ \ 2k / • \ 2fc_1 

% 7o Ve'/ 2 f-l y t J^8gx(i^) 1 ' 

Both sums and the integral in the above formula may be performed analytically leading to 

Hence we obtain 

Tito, 1) = -fclL - 3)£l(9 ) + *M {L _ 3 ) lo g2 + \?0 + ^i-^ (B.14) 



+0(iog 5 -*m+v)+---- 



The anomalous dimension for m > 2 is consequently given by 

7 (L,n) = f(g) (logn + | 7e + (*±i) - (L - 3) log 2) + 16 5 2 7l ( 5 , 1) 

4 5 -^2) iog^ + 6 5 (-l + log2) + (l-L) + ^- + ^ 
ir / g m — 1 2 

+ ^_^ + 0(l) (B15) 
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